Exact Relativistic Static Charged Perfect Fluid Disks 
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Using the well-known "displace, cut and reflect" method used to generate disks from given solu- 
tions of Einstein field equations, we construct static charged disks made of perfect fluid based on the 
Reissner-Nordstrom solution in isotropic coordinates. We also derive a simple stability condition 
for charged and non charged perfect fluid disks. As expected, we find that the presence of charge 
increases the regions of instability of the disks. 
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I. INTRODUCTION 

Axisymmetric solutions of Einstein field equations 
corresponding to disklike configurations of matter have 
been extensively studied. These solutions can be static 
or stationary and with or without radial pressure. So- 
lutions for static thin disks without radial pressure were 
first studied by Bonnor and Sackfield , and Morgan 
and Morgan j^] , and with radial pressure by Morgan and 



Morgan [3j. Other classes of static thin disk solutions 
have been obtained while stationary thin disks 

were studied in (3 _ [l3- Also thin disks with radial ten- 
sion 0| , magnetic fields and both electric and mag- 
netic fields [13 have been studied. An exact solution to 
the problem of a rigidly rotating disk of dust in terms of 
ultraelliptic functions was reported in |l4j . while mod- 
els of static relativistic counterrotating dust disks were 
recently presented in ^5[. Structures of nonaxisymmet- 
ric matter distributions of static charged dust were also 
recently studied The non-linear superposition of a 
disk and a black hole was first considered by Lemos and 
Letelier ^t|. Models of thin disks and thin disks with 
halos made of perfect fluids were considered in 0] . The 
generalization of the "displace, cut and reflect" method 
(Sec. |nj to construct thick static disks was considered 
by Gonzalez and Letelier [l^. For a recent review on 
relativistic accretion disks see pcj . 

Gravitationally bound stellar objects are likely to 
be positively charged due to the fact that the electrons 
are lighter than the protons so the formers can more eas- 
ily scape from the stellar object. One should expect to 
have an equilibrium situation when the electrostatic en- 
ergy of the electron (eV) is of the order of its thermal 
energy (kT) [2l|, then the scape would stop. Similar 
considerations can be found in Ref. |22j . 

The aim of this paper is to present a charged 
version of perfect fluid disks that were studied in |18j . 
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They are constructed with the well known "displace, cut 
and reflect" method applied to the Reissner-Nordstrom 
metric in isotropic coordinates. The stability of these 
new model of charged disks will also be considered. We 
present and extension of the Rayleigh criteria of stability 
[23j for the case of relativistic disks made of a charged 
perfect fluid. We find a simple condition to have gravi- 
tationally bounded disks. 

The paper is divided as follows. Sec. ITTI discusses 
the formalism which can be used to construct disks given 
a solution of Einstein-Maxwell equations, and the main 
physical variables of the disk. In Sec. IIIII Rayleigh in- 
spired criteria of stability criteria for the stability of 
charged and noncharged perfect fluid disks is established. 
The properties and stability of the charged disks are dis- 
cussed in Sec. IIVI Finally, in Sec. we discuss our 
results and make some considerations about the con- 
struction of charged perfect fluid disks with halos. The 
"displace, cut and reflect" method is reviewed in an Ap- 
pendix. 



II. EINSTEIN-MAXWELL EQUATIONS AND 
DISKS 

The isotropic metric representing a static spheri- 
cally symmetric space-time can be expressed as 



ds 2 = e^dt 2 - e A M(dr 2 + rW), 



(1) 



where dQ 2 = d6 2 + sin 2 8dip 2 . The same metric in cylin- 
drical coordinates (t, R, z, if) reads: 



by 



ds 2 = e^ R ^dt 2 - e A <*.*) (di? 2 + dz 2 + R 2 dip 2 ) . (2) 
The Einstein-Maxwell system of equations is given 
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where all symbols have their usual meaning. We use ge- 
ometric units with G = c = 1 . 

The method used to generate the metric of the 
disk and its material content is the well known "displace, 
cut and reflect" method (see the A ppe ndix) that was first 
used by Kuzmin j^] and Toomre |25l ] to construct New- 
tonian models of disks, and later extended to general rel- 
ativity (see, for example |6j,ll0|). Given a solution of the 
Einstein-Maxwell equation, this procedure is mathemat- 
ically equivalent to apply the transformation z — > \z\ + a, 
with a constant, on that solution. In the Einstein tensor 
we have first and second derivatives of z. Remember- 
ing that d z \z\ — 2§{z) — 1 and d zz \z\ — 2S(z), where ^{z) 
and S(z) are, respectively, the Heaviside function and the 
Dirac distribution, Einstein-Maxwell equations give us 



G^u = 87r(X^, m ' 



AirJ v 8{z), 



(7) 
(8) 



where T™ 11- is the electromagnetic tensor Eq. (0J, is 
the energy-momentum tensor on the plane z = and J v 
is the current density on the plane z = 0. For the metric 
PJl. the non-zero components of are 



Q\ = ^[-b zz +9 zz {bK + K + b%)\ 



(9) 



Qr = Qv = [- bZZ + 9 zz (b\ + &g + K)] , (10) 

where denote the jump of the first derivatives of the 
metric tensor on the plane z = 0, 



fytif 9^v,z\z=0+ 9^v,z\z=0~ i 



(11) 



and the other quantities are evaluated at z = + . The 
electromagnetic potential for an electric field is 



A M = (0,0,0,0). 



(12) 



Using Eq. (|12H and Eq. (JSJ , the only non-zero component 
of the current density on the plane z = is 



J' 



1 
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(13) 



where et M denote the jump of the first derivatives of the 
electromagnetic potential on the plane z = 0, 



(14) 



and the other quantities are evaluated at z = + . The 
"physical measure" of length in the direction d z for 
metric J5J is \/—g zz , then the invariant distribution is 
S(z)/y/—g zz . Thus the "true" surface energy density a 
and the azimuthal and radial pressures or tensions (P) 
are: 

<T = V=g^Qt, P = -V=9TzQ% = ~V z: gTzQ%. (15) 

Since J M = pU^, where C/ M = 8% ' / ' -Jgu, the "true" surface 
charge density p is 



P = \Z-g Z z9ttJ 



(16) 



III. STABILITY CONDITIONS FOR PERFECT 
FLUID DISKS 

One way to explain stability of static disks without 
radial pressure is to assume that the particles on the disk 
plane move under the action of their own gravitational 
field in such a way that as many particles move clockwise 
as counterclockwise (counterrotating model). With this 
assumption, stability of the matter on the disk can be 
associated with stability of circular orbits along geodesies 
(see [2(| for a detailed explanation). The usual stability 
criteria for circular orbits is adapted from the Rayleigh 
criteria of stability for a rotating fluid [23j . One finds that 
stability against small radial perturbations is achieved 
when 



u dh n 

h— > 0, 

ar 



(17) 



where h is the specific angular momentum of the circular 
orbit. 

In the case of perfect fluid disks the situation is 
somewhat different, since radial pressure can equilibrate 
the inward gravitational force and no counterrotating hy- 
pothesis is needed. Using T^ u '.„ = where T^ v is the sum 
of the energy-momentum tensor for a perfect fluid and 
Eq. ijljl. the equilibrium condition in the radial direction 
on the plane z = is given by 



-pe 



z u P 



.R- 



(18) 



The left side of Eq. I|18H can be interpreted as the gravi- 
tational force which equilibrates the pressure and electric 
forces that appear on the right side. Now suppose an ele- 
ment of fluid at radius R is displaced to R + AR keeeping 
P, a and p constant. The right side of Eq. ljT%|> becomes 

-e^ R+AR y 2 <l> >R (R + AR)p(R)-e^ R + A ^P R (R). (19) 

These "forces" should be compared with the right side of 
Eq. O at radius R + AR: 

- e v{R+AR)/2 ^^ R {R + AR)p(R + AR) 

-e^ R+AR) P R (R+AR). (20) 

To have stability expression (|19|l must be less then ex- 
pression (1201) ■ Expanding p(R + AR) and P,r(R + AR) 
around R, we get 



tR p,R + e v P RR < 0. 



(21) 



For an uncharged fluid, condition (|21|l reduces to P^rr < 
0. 

Note that this criterium of stability refers only to 
radial perturbations of the pressure and charge. It really 
gives us a condition to have gravitationally bounded sys- 
tems. In other words the disks do not explode, but they 
can collapse. In the general case, the study stability of 
gravitating systems reduces to the much harder problem 
of the study of the eigenvalue problem for a nontrivial 
elliptic operator |27j . 
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IV. CHARGED PERFECT FLUID DISKS 



tial associated to solution (|22() is 



We apply now the results of the previous sections 
to construct charged disks. The Reissner-Nordstrom solu- 
tion in Schwarzschild coordinates is given by 



dr 2 



1_ 2m + Q1 



(22) 



Au = ( -,0,0,0 



With the radial coordinate transformation 



r = r' I 1 



m + Q 
2? 



m — Q 
2? 



(23) 



(24) 



where m and Q are, respectively, the mass and charge of metric l(22() and Eq. (|23|l can be expressed in isotropic 
the black hole, and m > Q. The electromagnetic poten- coordinates (t,r',d,ip) as 



ds 2 = 



(m 2 -Q 2 ) 
IP? 



rdt 7 



(m+Q) 



(m-Q) 



(m + Q) 
2r' 



(m-Q) 



2r' 



(dr M + /W + r z sin' (9d</), (25) 



-4„ 



,0,0,0 . (26) 



Transforming Eq. (|25() - (|26|l to cylindrical coordinates, 
and using Eq. 1)15(1 and 1)16(1. we obtain a disk with sur- 
face energy density <r = a/m, equal radial and azimuthal 
pressures (or tensions) P = P/m and surface charge den- 
sity p = p/m where 



P = 



45 2 VP 2 + a 2 + 1-Q 2 
*■ [(1 + 2^-R 2 + a 2 ) 2 - Q 2 ] 2 ' 
25 1 - Q 2 



8Qa\/R 2 + ~a 2 



Wn 2 + a 2 {l + VP 2 + a 2 ) + 1-Q 2 



(27) 



w [(1 + 2V^ + 5 2 ) 2 - Q 2 ][l - Q 2 - 4(i? 2 + a 2 )} 

(28) 



(29) 



with i? = R/m, a = a/m and Q — Q/m. 



Eq. l(2"7|) shows that the disk's surface density is 
always positive (weak energy condition) for Q < 1. Pos- 
itive values (pression) for the stresses in azimuthal and 



radial directions are obtained if a > y 1 — Q 2 /2. The 

velocity of sound propagation V, defined as V 2 = is 
calculated using Eq. ((23 and Eq. (J2HJ: 



(1 - g 2 )[(l + 2y/^ 2 + a 2 ) 2 - Q 2 ][(l + 2^J R 2 + 5 2 ) 2 (1 - 4y/^ 2 + a 2 ) - Q 2 
[1 - Q 2 - 4{R 2 + 5 2 )] 2 [-3(l + 2yjR 2 + a 2 ) 2 + Q 2 (3 + 8 VR 2 + a 2 )} 



Fig-Hshows the curves of V 2 = 1 (solid curve) and 
of 5 = yl — Q 2 1 2 (dotted curve) where P changes sign 

as functions of the parameters a = a/m and Q = Q/m. 
Above the dotted curve, stresses are positive (pressure) 
for all R and above the solid curve, condition V 2 < 1 is 



also satisfied for all R. Thus, choosing values for a and 
Q that lie above the solid curve ensures that the entire 
disk will have pressures and sublumial sound velocities. 
We also note that with increasing charge the range of 
the cut parameter a that generates disks for which the 
conditions stated above are satisfied is enlarged. 



FIG. 1: Curve V 2 = 1 (solid curve) as a function of the 
parameters a = a/m and Q = Q /m for the charged disk of 

perfect fluid. Also a = u 1 — Q 2 /2 (dotted curve). 



Fig. El (a)-(d) shows, respectively, the surface en- 
ergy density a, pressures P, sound velocity V and charge 
density p with a = 1, Q = 0, 0.3, 0.6 and 0.9 as functions 
of R. As charge increases, the disks become less rela- 
tivistic for the same cut parameter; energy density and 
pressures are lowered and charge density becomes more 
concentrated near the disk center. 

Fig. Ola) is a graph of curves where Eg. I2T1 changes 
sign. The curves have been plotted only for ranges of 
parameter a where V 2 < 1 (Fig. . At the left of each 
curve, stability condition (|21[) is satisfied. Thus the disks 
are stable only in a small region near their centers. We 
also note that the charge decreases the radii of stability. 
The left side of Eq.|^is plotted in Fig. OJb) for the same 
parameters as in Fig. |21 



V. DISCUSSION 

We applied the "displace, cut and reflect" method 
on the Reissner-Nordstrom solution in isotropic coordi- 
nates and generated static charged disks made of perfect 
fluid. We also derived a simple criteria for the stabil- 
ity for charged and uncharged perfect fluid disks that 
is an extension of the Rayleigh criteria of stability for 
rotating fluids. The addition of charge decreases the en- 
ergy density and pressures near the disk's centers, while 
charge density in enhanced there. Furthermore, pres- 
ence of charge decreases the stable regions of the disks, 
that are reduced to small regions near the center even 
in the absence of charge. This is a rather different re- 
sult from that of our previous stability analysis of the 
uncharged perfect fluid disk 0] based on the stability 
study of circular orbits of test particles, where we found 
that the disks were completely stable for a ^ 1.016. It 
is worthwhile to note that the criteria used to study sta- 
bility of disks made of counterrotating matter is based 
on a particle consideration, whereas the stability criteria 



QTo" 




0.1 1 1 1 1 ' 1 

12 3 4 5 
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FIG. 2: (a) The surface energy density a, Eq. H27H . (b) the 
pressure P, Eq. (1281 . (c) the sound velocity V Eq. il->( )| and 
(d) the charge density p, Eq. 12911 for the disk with 5=1, 
Q — 0, 0.3, 0.6 and 0.9 as functions of R. 
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Q=0.3 
Q=0.6 
Q=0.9 




exact solutions is more involved than in canonical spher- 
ical coordinates. We believe that this search for exact 
solutions in isotropic coordinates is worthwile. 





FIG. 4: An illustration of the "displace, cut and reflect" 
method for the generation of disks. In A the spacetime with a 
singularity is displaced and cut by a plane (dotted line), in B 
the part with singularities is disregarded and the upper part 
is reflected on the plane. 
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APPENDIX 



S. 



FIG. 3: (a) Curves where Eq. 12111 defining stability changes 
sign as functions of radius R, cut parameter a and charge Q. 
Disks are stable at the left of each curve, (b) Curves of the 
left side of Eq. 12 H for the same parameters as in Fig. [5] 



derived for perfect fluid disks is based on collective phe- 
nomena. Therefore, the stability criteria derived in this 
paper seems to us to be more appropriate for the study 
of perfect fluid disks. 

In principle it is possible to add charged halos to 
the disks presented in this work by applying the "dis- 
place, cut and reflect" method to a static charged sphere 
of perfect fluid in isotropic coordinates (see the Ap- 
pendix). The halos could have charge of the same sign 
or opposite sign of the disk's charge such that the en- 
tire object could be neutral or have an arbitrary charge. 
Although there exist many exact solutions for charged 
fluid spheres in Schwarzschild coordinates (see, for ex- 
ample, |28j1. we have not found in the literature similar 
solutions in isotropic coordinates. This may be due to 
the fact that in these coordinates the Einstein-Maxwell 
equations for a charged perfect fluid are reduced to a sys- 
tem of highly nonlinear coupled second order differential 
equations for both metric functions v(r) and A(r) and 
for the electrostatic potential <j>{r), so the task of finding 



In this appendix we give an overview of the "dis- 
place, cut and reflect" method used to generate the met- 
ric and its material and electric content from a known 
solution of the Einstein-Maxwell field equations. The 
method can be divided in the following steps that are 
illustrated in Fig. 0] First, in a space wherein we have a 
compact source of gravitational field, we choose a surface 
(in our case, the plane z = 0) that divides the space in 
two pieces: one with no singularities or sources and the 
other with the sources. Then we disregard the part of 
the space with singularities and use the surface to make 
an inversion of the nonsingular part of the space. This 
results in a space with a singularity that is a delta func- 
tion with support on z = 0. The same procedure can 
be used with a static sphere of charged perfect fluid to 
generate charged disks with halos, as depicted in Fig. 
the sphere is displaced and cut by a distance a less 
then its radius. The part of the space that contains the 
center of the sphere is disregarded. After the inversion 
of the remaining space, one ends up with a charged disk 
surrounded by a cap of charged perfect fluid. If the inter- 
nal fluid solution is matched to the Reissner-Nordstrom 
metric Eq. I|25[l. the outer part of the disk will have the 
physical properties deduced in Sec. II VI while the proper- 
ties of the inner part will depend on the particular fluid 
solution. 
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FIG. 5: An illustration of the "displace, cut and reflect" 
method for the generation of charged disks with halos. In 
A the sphere of perfect fluid is displaced and cut by a plane 
(dotted line), in B the part that contains the center of the 
sphere is disregarded and the upper part is reflected on the 
plane. 
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